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We construct an effective low energy Lagrangian of gluodynamics which (i) satisfies 



all constraints imposed by the Renormalization Group; (ii) is scale and conformally 



Qh! invariant in the limit of vanishing vacuum energy density; (iii) matches onto the 



perturbative theory at short distances. This effective theory has a dual descrip- 
tion as classical gluodynamics on a curved conformal background. Color fields are 
dynamically confined, and the strong coupling freezes at distances larger than the 
glueball size. We also make specific predictions (in particular, on the N c dependence 
of glueball properties) which can be tested in lattice simulations of gluodynamics. 

I. INTRODUCTION 

A. Instability of the perturbative vacuum 

It has been known for a long time that the perturbative QCD vacuum is not the true 
vacuum of the theory. One way to see this is to examine the derivation of the asymptotic 
freedom jlj in the effective potential method. The real part of the one-loop potential of 



gluodynamics for a constant chromomagnetic field H reads 



ReV pert (F) = \H 2 + (9H) 2 ^- 2 (ln^f - ^) , (1) 



where \i is the renormalization scale, and b = lliV c /3. This potential has a minimum at 
H — H : © 

5fi? = A* e 6s2(Al) - (2) 

which is natural to identify with the perturbative vacuum. 

However it was soon realized that this perturbative vacuum is unstable. It is instructive 
to trace the origin of this instability in the effective potential method, which was pointed 
out in Ref. jjj. Consider the Landau levels of a particle of spin s and four-momentum p in a 



constant chromomagnetic field Hz 

p^p tl = 2gH{n + l/2)-2s z9 H, (3) 



where s z is a projection of the spin on the direction of the chromomagnetic field. The 
effective potential (JTJ can be calculated as [sj 

<jH 

I apz 

n=0 s„=±l 

Its real part yields (|T|l. while the imaginary part can be calculated as 



n TT „ OO 

Vp ert (iO = y d p z Y, E pgH{n + l/2-s z )+pl. (4) 

t: 7T J „_n „ LI 



OH f+9H . n^H 2 

lmV pert (H) = \ dp z ^p 2 z -gH-iO = (5) 

All 2 J-gH V 87T 

Therefore, the perturbative vacuum corresponding to the minimum of the perturbative po- 
tential (j2J is unstable. The instability is caused by the mode n = and s z = 1 (spin 
direction parallel to the field); note that n = corresponds to Landau level of the largest 



radius ~ 1/y/gH, i.e. to the infrarec 
QCD is ill-defined at large distances 
of other variables. 



5, 



region of the theory. This means that perturbative 
, and we may have to describe the theory in terms 



B. QCD in a cavity 



The breakdown of the perturbative approach (at least, at the one-loop level) has to 
happen at some critical value of the chromomagnetic field H c > Hq = An CD /g. This means 
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that weaker color fields cannot penetrate the physical vacuum, and the necessary condition 
for the applicability of the perturbative approach is that the energy density of the color field 
is sufficiently high: 



tX = ^~> 7^7— (6) 



?7T 327r 2 a s 

The condition (jHJ) means that the color fields can be properly defined only at distances 
smaller than -R con f ~ ^qcd- F° r perturbative theory to make sense, it has therefore be 
constrained within a cavity of radius -R con f, with appropriate boundary conditions. A possible 
realization of this idea is the MIT bag model [y] where the colored fields are required to vanish 
at the surface of a sphere. 

It is well known that the presence of boundary conditions leads to the emergence of 
Casimir vacuum energy ec ~ -R^tf ~ ^-qcd'i i n MIT bag model, it is represented by a "bag 
constant". It is also known that a theory in flat space-time in the presence of non-trivial 
boundary conditions can often be conveniently described as a theory in a curved background 
In this paper we will argue that such a description is possible for gluodynamics. We 
develop an effective theory which has the following dual descriptions: (i) classical gluody- 
namics in a curved conformal space-time background and (ii) gluodynamics in flat space-time 
coupled to scalar glueballs (which in this case play the role of dilatons saturating the cor- 
relation functions of the trace of the energy-momentum tensor). The representation of the 
effective theory in flat space-time appears quite similar to the non-topological soliton model 
of Friedberg and Lee jfj which describes quarks coupled to a scalar self-interacting field, and 
more generally to the approach outlined in Ref. we will return to the discussion of this 
topic later. 

It is clear that such an approach should have its limitations. Consider, for example, the 
dependence on the number of colors N c : the energy density of the gluon field e# ~ (AT? — 1), 
whereas the Casimir vacuum energy ec ~ ^qcd ~ N®. One therefore can expect that the 
effect of the boundary will diminish at large N c , and so the approach may not have a smooth 
N r — > oo limit. 



C. Renormalization group and low energy theorems 



The basic property of the perturbative effective potential is its invariance under the 
Renormalization Group (RG) transformations. We would like to preserve this fundamental 



property at all distances 
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j. For this purpose, we need to encode the properties of 
RG in a set of low energy theorems for the correlation functions of the trace of the energy 
momentum tensor. 

Let us sketch the derivation of these theorems, as they represent the guiding principle for 
the construction of our effective theory. Consider an expectation value of an operator O of 
canonical dimension d; it can be written down as 

s» id, 



(O) 



M e 



(7) 



On the other hand, the dependence of the QCD Lagrangian on the coupling is 



£ QCD = (-1/V)i>^i> , (8) 



13. 



14| | we can write down the 



/ dx <T{ 0(x) , F 2 (0) }) = - ^ (O). (9) 



where F = gF is the rescaled gluon field. Following Ref. 
expectation value of the operator O in the form of the functional integral and differentiate 
with respect to (— l/4g 2 (/x)) to get 

d_ 

Combining (J7J) and (jHJ) we obtain the relation [l3|, \l£ 

i lim / dxe i9X (0\T{O(x), ^ F 2 (0) }|0) connected = (O) (-4) . (10) 

This expression can be easily iterated by consequent differentiation like in Q to obtain a 
set of relations between Green's functions involving an arbitrary number of operators F 2 . 
We can rewrite those relations using expression for the scale anomaly in QCD in terms of 
the trace of the energy momentum tensor 9^ {d = 4) 

Substituting also 9fi for O we obtain the following set of low energy theorems for different 
Green's functions involving operator 0^(x): 

V 1 J ^ x . . . (OITI^H^O, • - • , ^:(^), ^(0)}|0) connected = (-4)*. (12) 

Eqs. f)l())l .(|12 jl show that although the scale symmetry of the classical Yang-Mills (JHJ) has 
been broken down by quantum fluctuations , there still remains a symmetry imposed by 
the invariance of the observables under the Renormalization Group. In the next section we 
are going to construct an effective Lagrangian which saturates (II 2j) at long distances and 
matches onto the one-loop perturbative effective Lagrangian at short distances. 
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II. EFFECTIVE LAGRANGIAN 



We start with the derivation of the effective Lagrangian using the mathematical trick 



suggested in llj. Consider the Yang-Mills theory on a curved conformally flat background 



in d dimensions. The background is given by the metric 

g^(x) = e^'^ (13) 

and the action by 

s = - jp f^xy/=ggr g ACT f; x P u a a . (i4) 

where g = detg M „. Recall that the classical Yang-Mills Lagrangian in flat spacetime is 
scale and conformally invariant only in four dimensions. On the contrary, it can be proven 



[llj that the theory on the curved background given by (|13|) .()14 |) is scale and conformally 
invariant in any number of dimensions d - this means that regularization does not bring into 
the theory (|14|) any dimensionful parameters. Upon regularization the action (|14|) acquires 
an additional term in d = 4: 



AS = l — I d 4 xe 2h 

Ag 2 



32 7T 



flUJ 



(15) 



The effective one-loop action of Yang-Mills field in the external constant conformally flat 
gravitational field is given by the sum of ()14j) and f|15j) : it is obviously scale and conformally 
invariant. The term (|15|) corresponds to the anomalous second term in the right hand 
side of (^). Therefore, the scale anomaly of QCD manifests itself in the theory defined 



by (|13|) and (|14j) through a term containing the axillary scalar field h ll| . without any 
dimensionful parameters. In a dual, and more conventional, flat space-time description the 
scale anomaly exhibits itself in the phenomenon of dimensional transmutation, which brings 
in a dimensionful parameter explicitly. 

The kinetic part for the field h(x) can be obtained in a manifestly scale and conformally 
invariant way using the Einstein-Hilbert Lagrangian for the one-loop effective Yang-Mills 
field 

S = Jd 4 xV^ R ~^ V g A<T *£x K, - e " h 6 ") > ( 16 ) 

where R is the Ricci scalar and G is some dimensionful constant; we substituted (|TT|) into the 
square brackets of JT3J). We can now use a well-known expression for the Riemann tensor 
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R^ u [25] to write down the dynamical terms for the field h(x) which obey the scale and 
conformal symmetry. Using (|T3*j) we get 

i?V=g = ^V=g = (17) 

Note that by writing (|T7jl we explicitly neglected terms of higher order in derivatives and 
constrained ourselves to the Einstein's gravity. This correspond to an expansion in powers 
of a slowly varying background field. 

The vacuum expectation value of the energy-momentum tensor reads 

<^> = -4|e v |. (18) 

The perturbative contribution to (jl8j) is given by Since the perturbative vacuum (J2J) 

is not stable, it is natural to assume that the dominant contribution to the energy density 
of the physical vacuum comes from non-perturbative modes. It is therefore convenient to 
separate the perturbative contribution to the 9£ in the following way: 

6* = expert.) -4|e v |; (19) 

(we will argue below (see (jUJ) that the physical vacuum is indeed independent of the 
value of the external chromomagnetic field.) Combining (|T6|) . (|T7j) and (|19|) we arrive at the 
expression for the effective one-loop action in the conformally flat gravitational field 



(20) 



where the new dimensionful constant m was introduced instead of G [lf| . 

At this point it is important to note that one can easily read the running coupling constant 
off (JT3J) and ()15|) . It can be seen that — e 2h plays a role of the familiar perturbative logarithm 
21n(g 2 //i 2 ). Hence our effective theory is applicable when q 2 < fi 2 . In the infrared region 
the perturbative expressions break down. However it is possible to remove the explicit 
dependence on the strong coupling from the effective action by performing the following 
redefinition in (|2()j) 

h — > h — 2 lm? , |e v | — > ^ 4 |e v | , m 2 — > q 2 m 2 , (21) 
where <^ 4 = g 2 6/32 tt 2 . Eq. (|21j) is just a change of mass unit. 
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Finally, we have to perform Legendre transformation of the action (}2T)j) to get the mini- 
mum of the effective potential at the minimum of the field \ which is canonically conjugated 
to the field h jllj]. The result reads 

C = ^ \ e* /2 (d»x? + |6 V | e* (1 - X) ~ \ (FP 2 +e*0-x)\ (F^)\ (22) 



where we included the factor y/— g in the definition of C The Lagrangian (}22|) defines our 
effective theory. At large distances (J22J) reduces to the low energy effective Lagrangian of 

Refs. HQ 

. Indeed, the Yang-Mills action is scale invariant in the external field h, which 
implies that ~ q 2 — > at long distances. We will investigate the region of applicability 



of the effective Lagrangian (j22j) in the following sections. 

The mathematical trick of putting the theory in curved space-time background which we 
used in derivation of f!22l) gives a simple way to keep track of all symmetries of the effective 
Lagrangian. However, we think it is also instructive to check how the effective perturbative 
potential (0) emerges from the Lagrangian f!22l) . The energy density 8 00 corresponding to 
(J22J) is given by 

e 00 (x) = M [{doX y + {diX f] e x/2 _ g oo | 6v | e x (1 _ x) 

+ (-F a0X F a \ + ig 00 (Ft) 2 ) (1 - e*(l - X )) , (23) 

where i — 1, 2, 3. Therefore the effective potential W in the constant chromomagnetic field 
H is 

W = jd?x !±H 2 - e*(l " X) (|^ 2 + k v |)} - (24) 

In strong chromomagnetic field H 2 ^> |e v | the energy density W reduces to the effective 
potential ((H). In this case x is n °t an independent degree of freedom, but rather a function 
of H. We calculate the corresponding momentum scale in the next section. The minimum 
of the functional W(H, x) is found from the following equations 

l- e x(l_ x ) = 0> X e x (\h 2 + |e v |) = 0. (25) 

Evidently, the minimum is at x = and the value of the W at the minimum is — |e v | 
independently of the value of the chromomagnetic field H. We conclude that the physical 
vacuum of the gluodynamics is described by one scalar field even in the presence of the 
applied chromomagnetic field. This justifies our assumption (|19)1. 
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It is seen from (I24|) that an increase of the color field H leads to the increase of the 
energy density of the system. Since H ~ g/r 2 (where r is the size of the system) the energy 
density decreases with r. At the same time the volume which the system occupies increases 
as r 3 . Therefore, we expect that there exists a static configuration with a finite size r$ such 
that the total energy of the system in minimal. This is analogous to the mechanism of bag 
formation in the Friedberg-Lee model However, the minimum of the effective potential 
in our model is independent of H and located at x = 0, while in the Friedberg-Lee model it 
depends on the density of the color sources. 

Note that we can read the one-loop behavior of the strong coupling right off the expression 
(£Q) for the effective potential. Indeed, the susceptibility of the vacuum in the strong external 
chromomagnetic field is (do not confuse fi(H) with the renormalization scale /z in ((H)) 

Recall that the beta function can be interpreted as a response of the system to the change 
of the external field. Namely, (j2*Uj) implies 

< 27 > 

From ()22j) it follows that jJ,{H) is independent of H at long distances, therefore (3 = 0. The 
strong coupling does not run if the effective theory is considered at the tree level. We will 
see in the next section that quantum corrections do not alter that conclusion. 

It remains to check that the vacuum at x = is stable. Let us recall that (£Q) is the real 
part of the perturbative effective potential. However the perturbative potential has also the 
imaginary part, as discussed above, which is due to the instability of the Landau level with 
n = and s z = 1, i.e. spin direction is parallel to the field. 

Let us now examine the properties of the Landau levels in our effective theory near the 
X = 0. The equation of motion of the dilaton field is 

^ d, {e x/2 d, X ) - M e */ 2 {d, X ? + X e* M + X #\ F%F*» = 0. (28) 
Expanding near the minimum we arrive at 

%X + m 2 (l + #Y) X = °- (29) 
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The corresponding Landau levels are 

= m 2 (l + 2^) • (30) 

It is seen that > for any H so that the instability does not develop in the effective 
theory we are discussing in this paper. 

Next, consider the trace of energy- momentum tensor which can be calculated directly 
from (f2"2*|) using 

W = (2^- - 5^C] + ^dle^\ (31) 

where the last term in the right hand side is the total derivative. Using equation of motion 
of the dilaton field ([28)1 one arrives at 

0£ = -4|e v |e* - X e x F« u F a ^. (32) 

By virtue of ()28j) one can clearly see that in the limit |e v | ^ the trace (|3*2*jl vanishes and 
the classical symmetries of the Yang-Mills Lagrangian are restored. One might be worried 
that in the expansion of (}3*2j) in powers of \ the term x^au F aiMU appears while it is absent 
in the Lagrangian (|22|). However it is easy to see that this term is canceled out by the pure 
dilaton contribution. Indeed, expanding the equation of motion (|2*5j) up to the quadratic 
terms in x one finds 

1 ^ pailv x = _ |ev| x + derivatiye + 0(x 2 } (33) 

It is important to check that 6£ satisfies the low energy theorems p2jl . We have noted 
just after (|2*2*j) that F" ~ q 2 . Hence it does not contribute to the 9% at q — > 0. The only 
remaining term in (|3*2*|) is the first one; it saturates the low energy theorems as was explicitly 



proved in |ll|. Thus, we have at our disposal the effective Lagrangian which (i) satisfies the 
constraints imposed by the Renormalization Group, (|32j) : (ii) the classical minimum x — 
of its effective potential (|2~%j) is the true stable vacuum of QCD; (iii) its scale and conformal 
symmetry is restored if |e v | = 0. 

III. QUANTUM FLUCTUATIONS AROUND THE PHYSICAL VACUUM 



Effective theory (|22|) is non-renormalizable. Let M be its UV cutoff (in the effective 
potential method, this is the scale which corresponds to the lowest Landau level). Quantum 
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fluctuations can develop only if there is enough kinematical space which is the case if m <C M. 
Let us define the perturbative expansion parameter A as 

We will see later in this section that perturbative series in powers of A is equivalent to 
the expansion of the Lagrangian (}2*2*|) in powers of Xi an d A indeed is the small expansion 
parameter in our effective theory. For the rest of this section we assume that A is small and 
prove this assumption in section ITVl 



A. Normalization of the energy-momentum tensor 

Let us first find the scale Mo at which our effective description breaks down; we will work 
in the leading order in A. The vacuum expectation value of the trace of energy-momentum 
tensor ()32j) is the physical observable and does not depend on a particular choice of degrees 
of freedom in the Lagrangian; its value is given by ()18j) . By virtue of (f3*2*|) it is equivalent to 
the requirement that 

4|e v |(l - e*> = (xe x F^ u F aflu ). (35) 

In the vacuum x — (jlSj) is obviously satisfied. Quantum fluctuations in general violate 
this requirement. However, since the effective Lagrangian (|22|) is formally divergent at short 
distances we have to impose an ultra-violet cutoff M . We will choose such a cutoff that 
(|18|) is satisfied. 

Expanding (J3*2j) to the order O(x ) we obtain a trivial result 

6jH = -4|e v | + Q(x). (36) 

At the order 0(x) Eq. (|33j) is satisfied due to At the next order 0(x 2 ) P°T> can be 

satisfied only for a particular choice of the cutoff Mq. Note that by the low energy theorems 
(fT2*j) (with n — 1) long distance contributions to expectation value of the operator 0£(x) can 
be expressed through the two point correlator S(g 2 ) defined as 

S(g 2 ) = i J d*xe^(0\T6Z(x)6Z(0)\0) = J da 2 ^ - , (37) 

where we have introduced the spectral density p(q 2 )- We find it more convenient to work 
with this correlator. The first reason is that the spectral density can be expressed in terms 
of physical states. The other one is that we know p(q 2 ) for the perturbative theory. 
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FIG. 1: (a) Pure dilaton contribution to the trace of energy-momentum tensor; (b) mixed dilaton- 
gluon contribution to the two-point correlator of the trace of energy-momentum tensor. 



To rewrite condition ()35|) in terms of two-point correlator we apply to it the low energy 
theorem (|12|) 

S(O) = -4<0£). (38) 

Thus, our requirement that (|35j) holds at the leading non-trivial order in A can be written 
as (see fl32l) 

H dil (0) + ~ mix (0) = , (39) 

where we separated the pure dilaton and mixed dilaton-gluon contributions. 
The pure dilaton contribution can be read from (|32j) : 

(^)dii = -4|e v |i( X 2 ) (40) 

which implies that (see Fig. ^a)) 

_ . . ., ,. . . m 2 1 f d 4 k i m 2 Ml , t , 

Sdii(O) = 8 |e v | < X 2 > = SKIt-t- / 777-777 = -77^ • ( 41 ) 



|e v | 2 J (2tt) 4 k 2 2tt 2 
Let us turn to the mixed gluon-dilaton contributions. The corresponding diagram is 
shown in Fig. db). Its imaginary part is 



1 '' % (2n) 2 5(k 2 1 )5(k 2 2 )(27Ty 5(k 1 + k 2 + qi + q 2 ) 



(2tt) 4 J (2tt) 4 

s / 2 \ 2 

o , 2 ^ / m 



(iV 2 - 1) f-7 , (42) 



140-48(2tt) 5V c 'Vie, 

where we neglected the mass of the dilaton m with respect to the cutoff M . S(g 2 ) can be 
calculated using dispersion relation with subtractions 

„ ( 2n f°° Pmix(q- 2 ) ^ 2 2 f°° Pmix(^ 2 ) ^Cr 2 



0~ 2 — Q 2 — i0 JO <T 4 
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FIG. 2: The leading order quantum correction to the gluon propagator at long distances. 



■i 



Pmix(o" 2 ) da 2 12 Z" 00 p mix (cr 2 ) da 



-1 / e ... - q , | j 



9 L ZuuZ^ ~2 77v\ " /„ ^2 ■ ( 43 ) 



r°° Pmix(^ 2 ) da 2 _ r°° p mix (v 2 ) da 

/„ a W ( a 2 _ q 2 _ z0 ) J 

The dispersion integral in the last line of (f4*HJ) is proportional to g 10 ln(— Mq + q 2 ). Conse- 
quently, 

S mix (0) = - T Pm ' x(0r 2 2) = - \ Pmix (M 2 ). (44) 
JO cr 4 



Formally, (|44|) gives the value of the non-vanishing subtraction constant in the dispersion 
relation. 

Substituting (j4~Tj) and (jjjj) into vacuum stability condition results in the equation 
determining the ultra-violet cutoff Mq of the effective theory 



2/3 

M 2 = 16 tt 105 1/3 (N 2 - l)- 1 / 3 ( ) . (45) 



B. Gluon polarization tensor 

We have argued that the vacuum expectation of the gluon condensate (fl"8|) is unchanged 
provided we had chosen the value of the cut-off according to (|4"Hj) . In that case the quantum 
corrections does not change the vacuum energy density which is completely saturated by 
classical solution. Now we would like to calculate quantum corrections to the strong coupling. 
To the leading order in A we have the tadpole diagram in Fig. El Introduce the scalar function 
II(g 2 ) as follows 

IV(g) = (<f<f - g 2 <ni%). (46) 
The tadpole diagram is given by 

- lw\C Wr^rro< ( - m ' q " ~ iry (47) 
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It can be calculated by performing Wick rotation and consequent integration over a four 
dimensional sphere of radius M . We neglect then the dilaton mass which gives contribution 
of higher order in A. The result of calculation of the diagram in Fig. [21 is 

n( " 2) = eJfci' 92£M «- (48) 

The quantum correction is constant. This means that the strong coupling freezes at long 
distances. 

The tadpole diagram Fig. El is leading order in A correction. However the higher order 
corrections deserve a special remark since they could be in principle logarithmically diver- 
gent at q 2 = Am 2 in which case those diagrams would dominate the polarization tensor at 
long distances. In the Appendix B we argue that all such logarithms appear in the product 
(q 2 — 4m 2 ) n ln(g 2 /4m 2 ), where n > and thus vanish at the end-point of the dilaton spec- 
trum. Also we check that the sub-leading diagrams are numerically small. Therefore, the 
conclusion of the previous section that the strong coupling freezes at long distances holds if 
such quantum corrections are included. 

Eq. (|35j) is the leading order contribution of vacuum fluctuations to the gluon polarization 
tensor. We can systematically develop the perturbation theory in A. The qualitative picture 
of the renormalization group flow can be obtained by simple dimensional analysis (see e.g. 
24j). Since the typical scale for mass is the cutoff M , the coefficients in front of the four 
terms in the rhs of (}22|) have the following behavior at different momentum scales p: (Mo/p) 2 , 
(Mq/p) 4 , (Mq/p)° and (Mo/p) 2 respectively. Thus, the only relevant term at low momenta 
is the second one, which is purely dilatonic term. This is a manifestation of the fact that 
the dynamics of the vacuum fields decouples from the colored sources. 

IV. MATCHING ONTO THE PERTURBATION THEORY 

One can express the strong coupling at the cutoff M as function of the parameters of 
the low energy Lagrangian. This can be achieved by matching the spectral density of the 
effective theory ()42|) with the spectral density of perturbation theory at Mq. In perturbative 
gluodynamics the anomalous trace of energy- momentum tensor is given by (lllj) . Then the 
calculation of the spectral density of the correlator (p?Tj) is straightforward |l7[ 

2 (ba s \ 2 2 N 2 - 1 f dSi f d% „ „ 2 
/W(S ) = I g^J 4 — ^— J — J ^(rffc - (ft ) 
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x(27r) 2 % 2 )% 2 2 ) (27r) 4 % + g 1 + g 2 ) 

q\ (49) 



'6a, V W - 1 ) i. 



\ 8tt / 2tt 

Since the spectral density is just the imaginary part of the correlator, it is clear that only 
mixed diagrams of ()44|1 contribute to the matching in the leading in A order (indeed, pan ~ 
m 4 ) 

Pmix(M 2 ) = p pcrt (M 2 ) . (50) 

Using Eq. (J43|) we obtain 

« S (M 2 ) = 1 . 6 f FA , -g 2 = Q 2 < M 2 , (51) 
6 JJV? - 1 



This equation shows that the small parameter of perturbation theory a s is matched onto 
the small parameter of our effective theory, A. 
At Q 2 > Mq the strong coupling runs as 



QCD 

where we introduced the familiar phenomenological constant constant Aqcd as 



4 TT 



A 2 =M 2 e^>. (53) 



Numerical estimations 



QCD sum rules analysis performed in 



13l |14| make it possible to estimate the non- 



perturbative scale inherent to the vacuum of gluodynamics, which appears quite hard: Mq = 
20 GeV 2 . Lattice calculations show that the lightest resonance in pure gluodynamics 
is the scalar glueball with mass m ~ 1.6 GeV. It is natural to identify this glueball with 
dilaton; it is interesting that this state appears to have a size much smaller than the sizes of 
glueballs with other quantum numbers 20j ■ In the approach followed in this paper this is 
a consequence of a large value of the cutoff scale Mq. From the vacuum stability condition 
(gHJ) we find |e v | ~ (0.58 GeV) 4 . By definition A = m/M ~ 0.36. Eq. (JSU) then implies that 
a s (M 2 ) ~ 0.33. The value of the A QCD follows from JE): A QCD ~ 0.79 GeV. 

In the world with light quarks the scalar glueball mixes with scalar qq meson Q]. The 
lightest scalar resonance is the cx-resonance which is a strong mixture of the glueball and 



15 

;he qq meson Q|. In this case the dilaton mass can be estimated as the mass of the a 
2]: m ~ 0.6 GeV. QCD sum rules give an estimate of the QCD vacuum energy density: 
|e v | = (0.24 GeV) 4 . From (|45)l we have Mq ~ 1.9 GeV. Other estimates can be done exactly 
as in the previous paragraph yielding A = 0.31, « s (Mq) ~ 0.35 and Aqcd = 0.26 GeV. To 
verify how good is this value from the phenomenological point of view we use ()52j) and find 

e a s (mz) — 0.12. This is in a reasonable agreement with the 



that at the Z-boson mass sea! 
data - see discussion in Ref. 



2l|. 



B. Dependence on N c 

Let us now discuss the dependence of our effective theory on the number of colors N c . 
Eq. (J55J) can be considered as an equation for the cutoff of the effective theory M as a 
function of the number of colors N c . Let us find Mq(N c ). It is convenient to introduce the 
dimensionless parameter a and have A rescaled as follows 



A = ^£ a = v^^T^^. (54) 
m v 8 m 



Then using and (J53J) Eq. (J53|) takes form 



a 



A = exp|--|. (55) 

Its solution is shown in Fig. E[a). We observe that the solution has two branches: one 
starting at the origin (0,0) and the second one starting at the point (0, 1). Both branches 
terminate at the critical point (a cr , A cr ) = (e -1 , e _1 ). To pick up the physical branch we note 
that by (p>4*|) a = when N c = 1. Thus, by ()51j) A = at this point. Therefore the physical 
branch is the lower one in Fig. Efa). In Fig. Efb) we represent it as a plot of the cutoff Mq 
versus the number of colors N c . The value of Aqcd — 0.8 chosen for this figure is such that 
a s (M|) = 0.12 at b = 11. The critical value of the parameter a corresponds to the critical 
value of A^ c . Using (|H1|) we find 

N*<Nr= ( Sp-e-iV + l, (56) 

(of course, the integer part of the right-hand-side must be taken) which yields = 3. 
When N c > N£ r our effective theory ceases to be valid. Indeed, M rapidly decreases with 
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FIG. 3: Numerical solution to if 55)). (a) Rescaled coupling A vs a; the critical point is (e _1 ,e _1 ). 
(b) Dependence of Mo on the number of colors in gluodynamics (Aqcd = 0.8 GeV); only the 
physical branch of the solution is shown. With a good accuracy Mq oc 1/N 2 . 

N c (approximately as 1/N 2 ) approaching the dilaton mass m. The values of the strong 
coupling a s (MQ r ) and the cutoff at the critical point are 

a s (M") = ^ , (57) 

and 

M? = . 8m ■ (58) 

We see that the effective theory breaks down at large N c , with the critical value 
(see ()56)) ). This can be readily interpreted if we recall that the perturbative vacuum energy 
density grows as |e v | pcrt ~ N 2 at large N c , whereas |e v | of the effective theory does not (see 
Fig. E^b) and ()45|l ). Thus, the effective theory breaks down at such large values of N c that 
the perturbative vacuum energy density cannot be matched onto the effective one. In the 
region N c < N£ r where we can use the effective theory, the iV c dependence of the value of 
the freezing strong coupling is given by a s {M^)b ~ X/N c ~ N c as can be seen from (|BT|) and 
Fig. lib). 

Since M decreases as N c increases the matching region is driven into the infrared where 
the perturbative expansion can no longer be trusted. Indeed, the dilaton spectral density 
vanishes at q 2 < Am 2 . On the contrary, the perturbative gluon spectral density is finite 
at arbitrary small but finite q 2 (see (|49Jl ). Although the dilaton effective theory takes into 
account the non-perturbative effects associated with the scale anomaly, it is not clear how 
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those effects are related to the color potential at long distances. The interplay between 
the dilaton low energy effective theory and the gluodynamics at large N c certainly deserves 
special study. 



V. CONCLUSIONS 

In this paper we constructed an effective low energy Lagrangian (|22j) of gluodynamics 
which (i) satisfies the constraints imposed by the Renormalization Group; (ii) its vacuum 
is stable; (iii) it is scale and conformally invariant in the limit of vanishing vacuum energy 
density |e v |; (iv) it matches onto the perturbative theory at short distances. Using this La- 
grangian we developed the perturbation theory of quantum fluctuations around the physical 
vacuum. Since the effective theory (|22|) is divergent when considered on a quantum level 
we must introduce an ultraviolet cutoff Mq. To calculate it we noted that classical config- 
uration of the dilaton field saturates the vacuum energy |e v |. Therefore the value of Mq is 
dictated by the requirement of vacuum stability - quantum fluctuations must not contribute 
to the vacuum energy density. This happens to be true only for a certain choice of M given 
by ([45)1 . In the kinematic region q 2 < Mq we developed a perturbation theory in a small 
parameter A = m/M Q and used it to calculate the leading (jjH| and next-to-leading (JB7j) 
order radiative corrections to the gluon propagator. We observed that the leading radiative 
correction to the gluon propagator is constant. We conclude that the strong coupling a s 
freezes at distances larger than the inverse cutoff 1/M n ; this behavior is consistent with the 

nn 

analysis of Refs. [21], 19] . 

By matching the spectral densities of the perturbation theory and of the effective one we 
determined the value of the strong coupling at the scale Mq in terms of the vacuum energy 
density |e v | and the glueball mass m, (fHTj) . Using QCD sum rules to estimate Mq ~ 20 GeV 2 
we calculate the Aqcd and then a s (mz); we found a reasonable agreement with experimental 
data. We consider this as an additional evidence that the typical scale of vacuum fluctuations 



of QCD is hard [l4, 



17 M 



2 

c ' 



We discussed the N c dependence of the theory. As N c increases M decreases as ~ 1/Af 
so that at some we have Mq < m and the quantum fluctuations of the dilaton field are 
no longer possible. The matching on the perturbation theory (JH0|) and (|51j) breaks down. 
Numerically, is found to be just above 3, so the effective theory ()22)1 is applicable to the 
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study the infrared behavior of SU (3) gluodynamics. 

One of our main results - the freezing of the strong coupling at long distances - has 
an elegant geometric interpretation. Recall that we derived the effective Lagrangian (|2*2*)l 
by for mally coupling Yang-Mills theory to the conformally flat gravity described by the 
field x (ill- This way the scale symmetry of Yang-Mills theory is restored at the cost of 
introducing a new field. At very short distances a s < 1 and the scale anomaly vanishes in 
usual perturbative gluodynamics. Effectively this means considering Yang-Mills theory in 
the flat space. At very long distances the theory resides in its physical vacuum \ — 0, see 
117)1 . which means that the spacetime is flat again. In between those extreme cases we can 
think of Yang- Mills field as a classical field propagating on a curved background. Indeed it 
has been found Ref. j3| that the coupling of the Yang-Mills theory on a curved background 
freezes at long distances. 

The physical picture which has emerged from our study thus corresponds to color fields 
dynamically confined within a cavity by the interaction with self-coupled scalar glueball 
fields. This interaction regularizes the theory in the infrared region, and leads to the freezing 
of strong coupling at large distances. 

It will be very interesting to study the properties of bound states in this "conformal bag 
model" . While we checked that the model does have the corresponding solutions, so far we 
have not succeeded in finding them analytically. 

A crucial test of the ideas presented in this paper can be performed on the lattice. Since 
r ~ 1/Mq corresponds to the size of the scalar glueball, and M Q decreases as a function of 
N c , we predict that the scalar glueball in 577(4) gauge theory will have a larger size than 
in SU(3). Unlike in 577(3) theory, where the scalar glueball was found to have the smallest 
size (see e.g. jsQ | ) , in SU(A) we expect all glueballs to have similar sizes. In contrast, in 
577(2) theory the size of the scalar glueball should become even smaller than in SU(3). 
These predictions can be tested directly by measuring the glueball formfactors (three-point 
correlation functions), or indirectly by measuring the two-point correlation functions of the 
scalar gluon operators and by checking at what distances they approach the perturbative 
behavior. If the lattice results in gluodynamics confirm the validity of the effective theory 
advocated in this paper, it will be worthwhile to include the light quarks by putting the 
classical QCD Lagrangian on the curved background. This could then substitute a consistent 
theoretical approach to the study of infrared behavior in QCD. 
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FIG. 4: Feynman rules for the dilaton effective theory. 
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APPENDIX A: FEYNMAN RULES FOR THE DILATON LAGRANGIAN 



In Appendix A we list the Feynman rules for the Lagrangian (|22jl up to the quadratic 
terms in \t see Fig- El Here a and b are the color indexes. We observe that dilaton graphs 
do not violate the color symmetry. This is seen of course directly from the Lagrangian (|22j). 
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APPENDIX B: HIGHER ORDER CORRECTIONS TO THE GLUON 

POLARIZATION TENSOR 

In Appendix B we argue that the higher order corrections to the gluon polarization tensor 
have no singularities at the end-point of the dilaton spectrum q 2 = 4m 2 . Let us consider 
the diagram Fig. El for example. 

We have 



m 



2 \ 2 1 r d*p d% d'k 2 UB / \ M-ig* 



i (p^q p - (qp)g 



e v \J 2\J (2vr) 4 (2vr) 4 (2vr) 4 v ^ y * ' p 2 



X i (pY - (qp)g^) 1 ^-— 2 ^-— 2 (2vr) 4 5(h + k 2 + p - q) . (Bl) 

Contracting Lorentz indexes and averaging over directions of p it can be shown that Il llL ,{q) 
has the same transverse structure as displayed in (|46j). Making contractions in the definition 
(|46|) we arrive at 

imn( ? ) = lining) = -L (f-Y -i- / d 4 fc ' f/V 



3g2 ^ 3g2 ^| £v |y 167r 7 7 ( 2 7t) ; 
x5(/c + p-g)5(p 2 )2(pg) 2 Wl 



4m 2 



M 2 



1 / m 2 \ 2 1 f „ , Q , -+, a . , / 4m 2 



3g 2 \|e v |/ 647T 3 



/ dM'lklulMjl-^, (B2) 



where M 2 is the external gluon virtuality. Denote t = M 2 . It is easily seen that 

t — M 2 

u p =\p\ = \k\ = t -^ r . (B3) 

Integral in ()B2jl over M 2 in the range 4m 2 < M 2 < t can be easily done giving somewhat 
lengthy result. Near the end-point of the spectrum the result of integration is 



2 



2 



(t - Am 2 f' 2 



Imn(t) « — ^ ; — rr^7- (B4) 

W Vl e vl/ 140 (67r) 3 mt 2 V 7 

The polarization tensor can be calculated using dispersion relation 

n (9 ) = - 4m")' - £ >2 (f _, 2 _ 80) ^ 4m y * ■ (B5) 

Dispersion relation can be applied only to a function which vanishes sufficiently fast at 
infinite radius in the complex plain of t. Therefore we apply it to a function ImlT(t) / (t— 4m 2 ) 5 
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FIG. 5: Next-to-leading order diagram contributing to the gluon polarization tensor. 

instead of ImII(£). This procedure corresponds to the subtractions 

4 1 

I% 2 ) -> II(g 2 ) - ^n«(4m 2 ) (g 2 - 4m 2 ) 1 . (B6) 



It follows from (|B4jl . (|B5j) and (JBJIJ) that 

11(g) oc (t - 4m 2 ) 9/2 -> 0, as t -> 4m 2 , 1^0. (B7) 



The term with I = is just the largest subtraction constant (cp. 

Therefore we can safely expand (|B2j) in powers of A. Integrating over M 2 and using 
dispersion relation ()B5jl with m = we obtain 

/^2x Q A ™ A i M 2 + Q 2 I 
nW) = W 1 " Q 2 24(4^ +C °° 3t - (B8) 

At Q 2 = Mq this contribution reaches its maximal value ~ A 2 and thus parametrically and 



numerically suppressed with respect to the leading result 

We can easily extend our argument to higher order diagrams. Indeed, the introduction 
of additional dilaton lines can only bring in a factors of M 2 / m 2 as can be seen from the 
gluon-dilaton vertex in the Appendix A and Fig. 1(b) of Ref. for dilaton self interactions. 
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